Observation of 7/;(4415) — > DD2(2460) decay using initial-state radiation 
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We report measurements of the exclusive cross section for e+e^ D^D^tv^ over the center- 
of-mass energy range 4.0 GeV to 5.0 GeV with initial-state radiation and the first observation of 
the decay i/)(4415) — > D^D^tt'^. From a study of the resonant substructure in 7/1(4415) decay we 
conclude that the '(/)(4415) D°D~-k+ decay is dominated by 7/;(4415) -> D!D2(2460). We obtain 
B(V'(4415) ^ J3°i3-7r+„__„J/B(V^(4415) ^ 00^2460) ^ D°D-7v+) < 0.22 at 90% C.L. The 
analysis is based on a data sample collected with the Belle detector with an integrated luminosity 
of 673 fb-^ 



PACS numbers: 13.66.Bc,13.87.Fh,14.40.Gx 

The V'(4415) resonance, the heaviest well established 
J^'-'' = charmonium state, was first observed thirty 
years ago by the MARK-I Q and DASP collabora- 
tions. Subsequently, although additional e~^e~ annihi- 
lation cross section measurements in the region of the 
-0(4415) were reported by the Crystal Ball [J| and BE- 
SII 0] groups, no update of its parameters were done 
until 2005, when a combined fit to Crystal Ball and BE- 
SII data was performed by Seth Q. Recently, the BES 
collaboration [5| reported new parameter values for the 
V'(3770), 0(4040), V'(4160) and 0(4415) resonances that 
are derived from a global fit to their cross section mea- 
surements. 

Despite the kinematic accessibility of ten open-charm 
strong decay modes for the -0(4415) its decays to any 
exclusive final states were not measured until last year, 
when a study of exclusive e+e" DD production via 
initial-state radiation (ISR) at Belle 0] revealed a small 
enhancement near the -0(4415) mass. A similar study @ 
of DD* production has no evident ?/'(4415) signal; the 
D*D* channel exhibits a small enhancement that may 
be attributable to the ?/;(4415), albeit at a slightly shifted 
mass value. 

In this Letter we report a measurement of the exclu- 
sive cross section for the process e+e" — > D°D~tt^ and 
the first observation of 0(4415) _DU2(2460) decay. 
It represents a continuation of our studies of the exclu- 
sive open charm production in the mass range where re- 
cently several new charmonium-like states were observed 
(y(4260) i, [3, F(4360), y(4660) X(4160) [13) 
decaying to either open- or closed-charm final states. Our 
study provides further information on the dynamics of 
charm quarks at these center of mass energies and on the 
properties of the -0(4415). The data sample corresponds 
to an integrated luminosity of 673 fb~^ collected with the 



Belle detector [l^l at the T(45') resonance and nearby 
continuum at the KEKB asymmetric-energy e'^e~ col- 
lider [13. 

We select e+e~ D'^ D'tt'^jisr signal candidates 
in which the D", D^ and 7r+ mesons arc fully recon- 
structed [ist - In general, the ^isr is not required to be 
detected and its presence in the event is inferred from a 
peak at zero in the spectrum of masses recoiling against 
the D^D^TT^ system. The recoil mass squared is de- 
fined as: M^^^{DDt:) = (E^c.m. - Edd-k)'^ - Pdd-k- Here 
^'c.m. is the initial e+e~ center-of-mass (cm.) energy, 
Edd-k and pddtt are the cm. energy and momentum 
of the D'^D~tt~^ system, respectively. To suppress back- 
ground two cases are considered: (1) the is out of 
detector acceptance in which case the polar angle for the 
D^D~Tr~^ system is required to satisfy | cos{6ijoD-Tr+)\ > 
0.9; (2) the fast is within the detector acceptance 
(I cos{6o<>D-7T+)\ < 0.9), in the latter case the jisr is re- 
quired to be detected with the mass of the D^D~Tr'^jjsR 
system greater than Ec.m. — 0.58GeV/c^. To suppress 
backgrounds from e+e^ DDiyk^jsr (n > 1) pro- 
cesses, we exclude events that contain additional charged 
tracks that are not used in D^ , D~ or tt"*" reconstruction. 

All charged tracks are required to originate from the 
vicinity of the interaction point (IP); we impose the re- 
quirements dr < 2 cm and < 4cm, where dr and \dz\ 
arc the impact parameters perpendicular to and along 
the beam direction with respect to the IP. Charged kaons 
are required to have a ratio of particle identification like- 
lihoods, Pk ~ ^K/i^K + ^-k) larger than 0.6. No 
identification requirements are applied for pion candi- 
dates. Kg candidates are reconstructed from tt^tt" pairs 
with an invariant mass within 10 MeV/c^ of the nominal 
Kg mass. The distance between the two pion tracks at 
the Kg vertex must be less than 1 cm, the transverse 
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flight distance from the interaction point is required to 
be greater than 0.1 cm, and the angle between the Kg 
momentum direction and the flight direction in the x — y 
plane should be smaller than 0.1 rad. Photons are recon- 
structed in the electromagnetic calorimeter as showers 
with energies greater than 50MeV that are not associ- 
ated with charged tracks. Pairs of photons are com- 
bined to form tt'^ candidates. If the mass of a 77 pair 
lies within 15MeV/c^ of the nominal 7r° mass, the pair 
is fit with a tt*' mass constraint and considered as a tt'^ 
candidate. D*^ candidates arc reconstructed using five 
decay modes: K^tt'^, K~K^, K^tt^tt^tt'^ , A'^tt+tt" 
and iC~7r+7r°. A ±15MeV/c^ mass window is used for 
all modes except for K^tt^tt^tt'^ , where a ±lGMeV/c^ 
requirement is applied ('^ 2.5a in each case). can- 
didates are reconstructed using the decay modes Kgir'^ 
and iir~7r+7r+. A ±15MeV/c^ mass window is used for 
all Z)+ modes. To improve the momentum resolution 
of D meson candidates, final tracks are fitted to a com- 
mon vertex with a mass constraint on the nominal D'^ 
or mass. Candidates in the D sideband region are 
selected for the background study and refitted to the cen- 
tral mass value of cither side of the sideband that has the 
same width as the signal window. To remove contribu- 
tions from e'^e" D~ D*'^^isr, -D^tt"*" combinations 
with invariant mass within ±10McV/c^ of the nominal 
mass are vetoed. 
The distribution of M'}^^{D°D-tt+) for the signal re- 
gion is shown in Fig. [1] (a). A clear peak correspond- 
ing to the process e+e^ D'^D^tt^^jsr is evident 
around zero. The shoulder at higher masses is due 

events, which include 
£)*o(-) _^ £)0(-)^o decays. To suppress the tail of such 
events we define the M^^j,(,(D°£'~7r+) signal region by a 
tight requirement ±0.7( GcV/c^)^ around zero. The po- 
lar angle distributions of D^D^-k^ and the mass spec- 
trum of the D^D~iT'^^isR combinations (when ^isr is 
detected) in the data shown in Fig. [1] (b) and Fig. [1] (c), 
respectively, are typical of ISR production and agree with 
the Monte Carlo (MC) simulation. The Mdod-^^ spec- 
trum after the application of all requirements is shown 
in Fig.[T](d). A clear peak is evident around the '0(4415) 
mass. 

The following sources of background are considered: 
(1) combinatorial background under the D^(D~) peak 
combined with a real D~{D^) coming from the sig- 
nal or other processes; (2) both and D~ candi- 
dates are combinatorial; (3) reflection from the processes 
e+e~ _D°D~7r+7rJJjjgg7/5ji'. with an extra Tr^j^igg in the 
final state, including -> D^^''> n^^^^^ decays; (4) 

reflection from the process e+e^ ^ D*'^ D^tt^^jsr, fol- 
lowed by D*'^ D'^^, with an extra soft 7 in the flnal 
state; (5) e+e~ D^D~tt^tt^ in which an energetic tt*' 
is misidentifled as a single Jisr- The contribution of 
background (1) is extracted using M^io and M^- side- 
bands that are four times as large as the signal region. 
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FIG. 1: The observed distributions of (a) M^^^{D'^ and 
(b) D'^D~tv'^ polar angles, (c) The D" D~ Tv'^'yisR mass spec- 
trum, (d) The resulting Mqoo-^+ spectrum. Histograms 
show the normalized contributions from the M^o and 
sidebands. The selected signal windows are indicated by ver- 
tical lines. 



These sidebands are shifted by BOMeV/c^ (20MeV/c2 
for the — > K^tt^tt^tt^ mode) from the signal re- 
gion to avoid signal over-subtraction. Background (2) is 
present in both the M^o and M^- sidebands and is, thus, 
subtracted twice. To account for this over-subtraction 
we use a 2-dimensional sideband region, where events 
are selected from both the M^o and the Mjj- side- 
bands. Backgrounds (1-2), ~ 15% of the signal, are 
subtracted from the signal-region £'°Z)~7r+ mass spec- 
trum. The dominant component of backgrounds (3-4) is 
suppressed by the tight requirement on M^^^{D^ D~ n^) . 
The remaining background from (3) is estimated by 
applying a similar full reconstruction method to the 
isospin-conjugated processes e^e^ D'^ D'^'n'^T^^issHSR 
and e+e~ — + D^D'tt^tt^^^^jisr- Here the absence 
of additional charged tracks in the event is not re- 
quired. Because of the charge imbalance in the D'^D'^tt'^ 
and D'^D^TT^ final states, only events with a missing 
extra -k^^^s can contribute to the M'^g^{D^ D'^ir^) and 
M^^^{D^ TT^) signal windows. To extract the level 
of background (3), the D°D°7r+ and D+Z)^7r+ mass 
spectra are rescaled according to the ratio oi and 
Z?" reconstruction efficiencies and a factor of 1/2 due to 
isospin. The averaged rescaled D^D^n^ and D^D^tt^ 
mass spectrum, which is small (~ 2% of the signal), is 
subtracted from the signal-region AIdod-^+ spectrum. 
Background (4) is also estimated from the data assum- 
ing isospin symmetry. We measure the process e~^e~ — > 
D*^ D^n+^iSR {D*~ D°7r~) applying a similar full 
reconstruction method. The shape and normalization 
of this spectrum after efficiency correction is then used 
to generate a MC sample of e+e^ D*^D^Tr^^isR 
(D*° — > -D°7) events. From this MC sample, tuned to the 
data, the contribution from background (4) is estimated 
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to be less than ^ 1.5 events/bin, which is subtracted 
from the signal Mdojj-^+ spectrum. The contribution 
from background (5) is determined from reconstructed 
e+e" Z)°Z)~7r+7r° events in the data found to be neg- 
ligibly small and taken into account in the systematic 
error. 

The e+e" — > D^D~'k^ cross section is extracted from 
the background-subtracted D°D~Tr~^ mass distribution 
following the procedure described in [1], taking into ac- 
count the differential ISR luminosity and the total ef- 
ficiency, which is found to have linear dependence on 
Md°d--k+ ■ The resulting e"'"e^ D'^D~tt^ exclusive 
cross section is shown in Fig. [2l Since the bin width is 
much larger than resolution, no correction for resolution 
is applied. 
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FIG. 2: The exclusive cross section for e^e — > Z)"Z) ■k'^ . 
The dotted line corresponds to the nominal mass of the 
V'(4415) [13]. 

The systematic errors for the a{e^e~ Z?°_D~7r+) 
measurements are summarized in Table U The system- 



TABLE L Contributions to the systematic error on the cross 
sections, [%]. 



Source 


D''D-n+ 


Background subtraction 


±5 


Cross section calculation 


±5 


B{D) 


±4 


Reconstruction 


±6 


Kaon identification 


±2 


Total 


±10 



atic errors associated with the background (1-2) subtrac- 
tion are estimated to be 2% due to an uncertainty in the 
scaling factors for the sideband subtractions. It is esti- 
mated using fits to the M]jo and Md- distributions in 
the data with different signal and background parame- 
terization. Backgrounds (3)-(5) are also subtracted us- 



ing the data and only the uncertainty in the scaling fac- 
tors for the subtracted distribution is taken into account. 
Uncertainties in backgrounds (3)-(5) are estimated con- 
servatively to be smaller than 4.5% of the signal. The 
systematic error ascribed to the cross section calculation 
includes a 1.5% error on the differential luminosity and 
4.5% error in the total efficiency fit. Another source of 
systematic error comes from uncertainties in track and 
photon reconstruction efficiencies (1% per track and 1.5% 
per photon). Other contributions come from the uncer- 
tainty in the kaon identification efficiency and the abso- 
lute and D~ branching fractions [l7j . 

To study the resonant structure in ?/)(4415) decays, we 
select D^D~iT^ combinations from a ±100MeV/c^ mass 
window around the nominal ?/'(4415) mass (toi/,(4415) = 
4.421 GeV/c2 [l^). A scatter plot of M{D~ii+) vs. 
M(Z)°7r+) and its projections onto both axes are shown 
in Figs. [3] (a),(b) and (c), respectively. Clear signals 
for the 75^(2460)° and 15^(2460)+ mesons are visible 
in these plots. We expect positive interference between 
the neutral Z)"!) J (2460)" and the charged D- Dl{2Am)+ 
decay amplitudes leading to the same D^D^-k^ final 
state for the decay of C = —1 state, and the scat- 
ter plot evidently agrees with this expectation. Be- 
cause of the interference we do not study D°Z32(2460)° 
and Z3~ £'2(2460)"'" final states separately and define 
the signal interval for the _DD2(2460) combinations 
as \Md-^+ - "i"D5(2460)nl < 50 McV/c^ or |M£,o^+ - 
™D5(2460)+I < 50MeV/c2 {m-^.^^^Q^y, = 2.461 GeV/c^, 
(2460)+ = 2.459 GeV/c2 [13). 

We perform a separate study of e'*"e~ D£'2(2460) 
and e+e" ■D(£'7r)„^„;o.(2460) ■ The Mood-^+ spec- 
trum for the £'£'2(2460) signal interval is shown in 
Fig. [IJa). A clear peak corresponding to ■i/;(4415) — > 
£'!D^(2460) decay is evident near the £'75^(2460) thresh- 
old. Assuming = 2+ for !D2(2460) mesons [13], the 
signal of -0(4415) £'U5(2460) should be described 
by a d-wave relativistic Breit-Wigner (RBW) function. 
However, to compare mass and width of the obtained 
-0(4415) signal with the corresponding -0(4415) reso- 
nance parameters measured in the inclusive study [l7j . 
we perform a likelihood fit to Mdod-^+ distribution 
with the £'£'2(2460) signal parameterized by an s-wave 
RBW function. The fit with the d-wave RBW function, 
which would be appropriate for the studied decay chan- 
nel, is found to be unstable with respect to variations in 
the background parameterization. To account for back- 
ground and a possible non-resonant D^D~-k~^ contribu- 
tion we use a threshold function y^A/ — m^) — "ni-Q* ^2A&a) 
with a floating normalization. Finally, the sum of the 
signal and background functions is multiplied by the 
mass-dependent linear efficiency function and differen- 
tial ISR luminosity. The fit, shown as a solid curve 
in Fig. HKa), yields 109 ± 25(stat.) signal events. The 
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FIG. 3: (a) The scatter plot of D"7r+ vs. D°-k+ for the signal 
region in the data for |M£,0£,-^+ — m^(44i5)| < lOOMeV/c^. 
{h)D~-K'^ and (c) D^tt'^ mass projections. Histograms show 
the normalized contributions from the Mjjo and M^- side- 
bands. The selected mass windows are illustrated by vertical 
and horizontal lines. 



significance for the signal is obtained from the quan- 
tity — 2 ln(£o/'Cmax), where £max is the maximum like- 
lihood returned by the fit, and Co is the likelihood 
with the amplitude of the Breit-Wigner function set to 
zero. Taking the reduction in the number of degrees of 
freedom into account, we determine the significance of 
the ^(4415) signal to be ~ 10a. The obtained peak 
mass m^(44i5) = (4.411 ± 0.007(stat.)) GeV/c^ and to- 
tal width Ttot = (77 ± 20(stat.)) MeV/c^ are in good 
agreement with the PDG [l7[ values, the recent BES re- 
sults and predictions of Ref. (3]. The peak cross 
section for e+e" ?/'(4415) ^ ZJ^Dj (2460) process at 
Ec.m. ~ ™?/;(44i5) is Calculated from the amplitude of the 
RBW function in the fit to be cr(e+e~ -0(4415)) x 
S(V'(4415) ^ DDi,{2A60)) x ^(75^(2460) ^ Dn+) = 
(0.74±0.17±0.08)nb. Here the systematic uncertainty is 
obtained by varying the fit range, histogram bin, param- 
eterization of the background function and efficiency. Us- 
ing cr(e+e- -> V(4415)) = 127r/TO| ^^^^^^ x (Tee/Ttot) we 

calculate the S(V'(4415) ^ DD^{24,60)) x S(i5*(2460) ^ 
Dtt+) = (10.5 ± 2.4 ± 3.8)% using the V(4415) parame- 
ters from the PDG [13 and (19.5 ± 4.5 ± 9.2)% for the 
■!/'(4415) parameters from Ref. 

The shape of the Mdod-^+ spectrum with the 
£'£'2(2460) signal excluded, shown solid curve in 
Fig. m^b). is relatively featureless and, in the 0(4415) 
mass window, is consistent with the combinatorial back- 
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FIG. 4: (a) The M£,od-^+ spectrum for the 00^(2460) sig- 
nal region. The solid curve represents the result of the fit 
described in the text. The threshold function is shown by 
the dashed curve, (b) The Mjyoi}-^+ spectrum outside the 
DDI (2460) signal region. The solid curve is a fit with a 
second-order polynomial background function. The dashed 
curve shows the upper limit on the ■0(4415) yield at 90% C.L. 
Histograms show the normalized contributions from Adjjo and 
M]j- sidebands. 



ground. The fit to this distribution with the -0(4415) 
signal function (an s-wavc RBW function with mass and 
width fixed to the result of the previous fit) and a second- 
order polynomial function yields a negative ^(4415) sig- 
nal that is consistent with zero. We calculate the upper 
limit on the 0'(4415) yield to be 18 events at 90% C.L. We 
assume a phase-space-like -0(4415) D^D^tt^ decay 
and calculate an upper limit on the ratio of the branch- 
ing fractions of -0(4415) decays to non-resonant D^D^tt^ 
and (2460). This assumption is not exactly correct 
since at least one meson pair has to be in a relative p- 
wave state. However, it provides a conservative estimate 
of the efficiency of the 13 1)2(2460) veto. In this way we 



obtain S(V'(4415) D°D-n, 



non — resonant 



)/S(0(4415) 



£>i:»^(2460) D°D-n+) < 0.22 at 90% C.L. 

In summary, we report first measurements of the 
e+e^ ^ D'^D^tt^ exclusive cross section over the center- 
of-mass energy range from 4.0 GeV to 5.0 GeV with ini- 
tial state radiation. Aside from a prominent 0(4415) 
peak, the cm. energy dependence of the e+e^ 
D^D~t:'^ cross section has no evident structure. From 
a study of the resonant structure in 0(4415) decay 
we conclude that the -0(4415) D^D^n^ process is 
dominated by 0;(4415) DDli'iA&Q). The mass and 
width of the 0(4415) state arc found to be (4.411 ± 
0.007(stat.)) GeV/c2 and (77±20(stat.)) McV/c^, respec- 
tively. 
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